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Abstract 

We study a mathematical model describing the dynamics of a pluripotent stem cell 
population involved in the blood production process in the bone marrow. This model 
is a diflerential equation with a time delay. The delay describes the cell cycle duration 
and is uniformly distributed on an interval. We obtain stability conditions independent 
of the delay. We also show that the distributed delay can destabilize the entire system. 
In particularly, it is shown that Hopf bifurcations can occur. 

Keywords: Blood production system, stem cells, delay differential equations, stability, Hopf 
bifurcation. 

1 Introduction 

Blood production process, called hematopoiesis, is one of the major biological phenomena 
occurring in human body. It takes place in the bone marrow where pluripotent stem cells 
give birth to mature cells. After ejecting their nuclei, these cells enter the bloodstream and 
become blood cells. 

According to the study of Burns and Tannock 4J , the population of pluripotent stem cells 
can be divided into two distinct groups: quiescent cells and proliferating cells. Mathematical 
models describing the dynamics of this cell population have been studied since the end of the 
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seventies, in particularly by Mackey [S1[TU]. We refer to the review articles by Haurie et al. 
[5] and Mackey et al. 12J for further study and more references on this topic. More recently, 
Pujo-Menjouet et al. pjj and Pujo-Menjouet and Mackey [T^j proved the existence of a 
Hopf bifurcation for the hematopoiesis model proposed in [H]. In all these works, the authors 
assumed that the proliferating phase duration is constant. Mathematically, this means that 
the delay in their models is a discrete delay. However, experimental data (see Bradford et al. 
[3j) indicate that cells do not spend the same time in the proliferating phase. 

In this paper, taking into account this assumption, we assume that the delay (or pro- 
liferating phase duration) is uniformly distributed on an interval. The main objective is to 
investigate the effect of time delay on the dynamical solutions. It is shown that there exist 
some critical values of time delay such that a local Hopf bifurcation occurs at the non-trivial 
equilibrium. 

The paper is organized as follows. In section [21 we present our model, which is given in 
equation ([T]) . In section O we derive stability conditions for the two equilibria of equation 
([!]) which do not depend on the delay. We show the existence of Hopf bifurcations at the 
non-trivial equilibrium in section [H A brief discussion is given in section [S] 



2 The Model 

Pluripotent stem cells can be either in a resting phase, also known as Go-phase, or in a 
proliferating phase. In the resting phase, they can die at a constant rate 6 > 0, which also 
includes the cellular differentiation, or be introduced in the proliferating phase at a rate /?. 
According to the work of Sachs [T^ , (3 is assumed to depend on the resting phase population. 

In the proliferating phase, which is in fact the so-called cell cycle, pluripotent stem cells 
are committed to divide and give birth to two daughter cells at the end of this phase. The 
two daughter cells enter directly the resting phase and complete the cycle. We assume 
that proliferating cells divide according to a uniform law / on an interval [Tmin,T] with 
< Tmin < T < +00. This assumption comes from the fact that, even if only a little is 
known about phenomena involved in hematopoiesis, there are strong evidences (see Bradford 
et al. indicating that cells do not divide at the same age. The function / is then defined 
by 

if r e [t„ 



fir) 



0, otherwise. 



Let x{t) denote the pluripotent stem cell population density (cells/kg) at time t >Q. It 
satisfies the nonlinear delay differential equation 

x'{t)^-{5 + l3{x{t)))x{t) + ^^ / (}{x{t ~ r))x{t - r)dr. (1) 

The first term in the right-hand side of equation ([T]) accounts for the cellular loss due to 
mortality and cellular differentiation, 5x{t)^ and introduction in the cell cycle, (3{x{t))x{t) . 
The second term is for the division of proliferating cells into two daughter cells during mitosis. 
Proliferating cells are in fact resting cells introduced in the proliferating phase one generation 
earlier, so that the quantity j3{x{t — r))x{t — r) appears with a time delay. The factor 2 is, 
of course, for the division of each proliferating cell into two daughter cells. 

In the following, the rate of reintroduction in the proliferating compartment /? is taken 
to be a monotone and decreasing Hill function, given by 

nn 

p(x) = Po- for a; > 0. 
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The coefficient /?o > is the maximum rate of reintroduction, > is the Go-phase popula- 
tion density for which the rate of re-entry /3 attains its maximum rate of change with respect 
to the resting phase population, and n > describes the sensitivity of (3 with changes in the 
population. This function was firstly used in hematopoiesis models by Mackey [1] in 1978. 

In [9] and [llj . Mackey gave values of the above parameters for a normal human body 
production. These values are 

(S = 0.05d"\ /3o = 1.77d"^ and n = 3. (2) 

The value of 9 is usually 9 = 1.62 x 10* cells/kg. However, since we shall study the qualitative 
behavior of the pluripotent stem cells population, the value of 9 is not really important and 
could be normalized without loss of generality. 

Now if we consider an initial continuous nonnegative function ip defined on [— r, 0], then 
the equation ((l|) has a unique continuous and nonnegative solution x'''(t), defined for t > — r, 
such that 

x'^(s) = (^(s) forse[-T,0]. 

This can be obtained by using the results in Hale and Verduyn Lunel [?]■ 

Notice that equation ([T]) has at most two equilibria, the trivial equilibrium a; = and a 
non-trivial positive equilibrium x = x* . The trivial equilibrium always exists and corresponds 
to the extinction of the population. 

Proposition 2.1. Equation ([Ip has a non-trivial positive equilibrium x = x* if and only if 

/3o>S>0. (3) 

In this case, x* is explicitly given by 




Proof. Let x* be an equilibrium of equation ([T]). Then x* satisfies 

x*{l3{x*)-S) =0. 

Consequently, equation ([1]) has a non-trivial equilibrium if and only if the equation 

/3(x*) 

has a non-trivial solution. Since the function /3 is decreasing and positive with /3(0) = /3o, 
then equation ([1]) has a non-trivial equilibrium if and only if condition Q holds. □ 

In the next section, we shall study the stability of the two equilibria of equation P]). 



3 Stability 

Throughout this section, we are interested in the stability of the equilibria of equation ((T|) , 
in particularly the stability of the non-trivial equilibrium x = x* . We start by giving a result 
on the global stability of the trivial equilibrium of ([l]). 

Theorem 3.1. The trivial equilibrium x = of equation (Qp is globally stable if 

Po<5. 
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Proof. The proof uses a similar technique employed by Adimy and Crauste [1 . It is based 
on the construction of a Lyapunov functional. 

Denote by C"*" the space of all continuous nonnegative functions on [— r, 0]. Let B be the 
function defined by 



B{x) = / /3(s)s ds for a; > 0. 
Jo 

Consider the mapping J : C+ [0, +oo) defined, for ip e C+, by 

= Bi^rn + — i / r ((3{p{a)Ma)) \adr. 

Then, 

= ^(0)/3(^(0))^(0) + ^— r (/3(^(0))^(0))' - (/3(^(-r))^(-r))'dr. 



Since 

^(0) = -(<5 + /3(^(0)))^(0) + — ^ / (3{^{-r)M-r)dr, 

we obtain that 

= -{5 + 0{ip{O)))p{cp{O))cp^O) + ^ r (/3(^(0)M0))'dr 

— r (/3(¥>(0)MO)-/3(^(-r)M-r))'dr. 

Hence 

j(^)<-(5-/3(^(0)))/3M0))v:>(0)2. 
Let a be the fmiction defined, for a; > 0, by 

a{x) = [5 - I3{x))l3{x)x^ . 

Assume that /3o < 5. Since /3 is a decreasing function, it follows that the function x ^ 6 — (3{x) 
is positive for a; > 0. Hence, a is nonnegative on [0, +oo) and a{x) = if and only if a:: = 0. 
Consequently, the mapping J is a Lyapunov functional when /3o < We then deduce that 
the trivial equilibrium of ([1]) is globally stable. □ 



The result in Theorem 13.11 describes the fact that when a; = is the only equilibrium of 
p^, the population is doomed to extinction except when (3q — S. 

Now we focus on the stability of the positive equilibrium a; = a:* of equation ([1]). To 
ensure the existence of the equilibrium x = x* , we assume that condition ^ holds; that is, 

Po>S>0. 

We do not expect to obtain conditions for the global stability of a; = a;*. However, local 
stability results can be obtained by linearizing equation ([T]) about a;*. Set 

The linearization of equation ([1]) at x* is 

2/3* 



X 



'(t) = -{S + I3*)x{t) + '- / x{t - r)dr. 



4 



M. Adimy, F. Crauste and S. Ruan 



A model of pluripotent stem cell dynamics 



The characteristic equation of ((T|) is given by 



A(A) —X + S + P* / e'^^'dr = 0. (5) 



We now state and prove our first resuh on the stabiHty of a; = cc*. 
Theorem 3.2. Assume that 

n—^ < 1. 

Then the non-trivial equilibrium x = x* of equation ([7]j is locally asymptotically stable. 

Proof. We show, in fact, that x* is stable when /3* > 0. By the definition of f3* given by ([4]), 
it follows that ^ 

/3* > if and only if "-"5 — — 1- 

Po 

So we assume that (3* > 0. 

We first assume that A(A), given by (O, is a real function. Then, A(A) is continuously 
differentiable and its first derivative is given by 

(A) = 1 + r re-^^dr. (6) 



One can see that dA/dX is positive for A G M as soon as (3* > 0. Moreover, 
lim A(A) = — cxD and lim A(A) = +oo. 

A — > — oo X — ' + 00 

Consequently, A (A) has a unique real root Aq. Since 

" 1% 



A(0) = ,5-/3* = - — ) > 0, 



we deduce that Ao < 0. 

Now, we show that if A is a characteristic root of equation then Re{X) < Aq. By 
contradiction, we assume that there exists a characteristic root X — ^ + iuj of equation ([5|) 
such that /i > Aq. By considering the real part of A(A), we obtain that 



li + S + f3* / e-''' cos{ujr)dr ^0. 



M - Ao = / (e-'^'^ cos{u;r) - e~^°'')dr < 0. 



Consequently, 



This yields a contradiction. We conclude that every characteristic root A of is such that 
Re{X) < Aq. Hence, all characteristic roots of ([5]) have negative real parts and the equilibrium 
X = X* is locally asymptotically stable. □ 

When /3* < 0, that is, when 

Kn——, 

Po 

the stability cannot occur for all values of Tmin and r. In particularly, we shall show that a 
Hopf bifurcation can occur (see Theorem l4.ip . However, we can still have the stability of the 
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non-trivial equilibrium x = x* for values of n, [3q and 5 if n(/3o — S)//3o is not too large. This 
will be considered in the next theorem. 

To present the results, without loss of generality we assume that 

We want to point out that the results we are going to show remain true when Tmin > 0, but 
the proof is more complicated. 

Define a function if, for cc > 0, by 

Kix) = ^ (7) 
and let xi be the unique solution of the equation 

xi = tan(a;i), a;ie(7r, — ). 

Set 

Mo '■— cos(a;i) £ (—1, 0). 

Then K'{xi) ^0 and 

Mo = K{xi) = inmK{x). (8) 

a;>0 

We have the following local stability theorem. 
Theorem 3.3. Assume that 

/3o 1 - 2mo ^ ^ 

Then the non-trivial equilibrium x = x* of equation ^ is locally asymptotically stable. 
Proof. Let us assume that © holds. Then /3* < 0, (5 + /3* > and 

6 + 13* 



2(3* 



< Mo. (10) 



By contradiction, assume that there exists a characteristic root X = fi + iuj oi ^ with /i > 0. 
Then, 



23* /"""^ r 
/Li = -((5 + + — — / e oj cos(r)dr 

Jo 



and 



2P* 



w = / e sin{r)dr. 

Jo 

Integrating by parts, we obtain that 

2/i = -{S + (3*) + 2l3*e-f^K{ujT). 

Consequently, 

^l<-{5 + 13*) + 2(3*e-^'''K{ujT). 

If LOT is such that 

sin(a;T) > 0, 
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then 

2f3*K{ujT) < 

and 

fi < ~i6 + f3*) < 0. 

So we obtain a contradiction. 
Similarly, if lut is such that 

sin(a;T) < 0, 

then, from ^ and (fTU]) . we deduce that 

— — < K(ijJT). 

2(3* - ^ ' 

It implies that 

6 + 13* > 2I3*K{lut) > 2/3*e~^^if(wT). 

Therefore, 

H<-{d + (3*) + 2I3*K{ujt) < 0. 

Again we obtain a contradiction. Hence, all characteristic roots A of (O are such that 
i?e(A) < 0. 

Now, we assume that ([5]) has a purely imaginary characteristic root A = iw. Then lo and 
r satisfy 

Kiivr) = ^±^. (11) 

Using ([S]) and PH)) . we obtain a contradiction. Consequently, pT|) has no solution and 
equation (O does not have purely imaginary roots. We conclude that all characteristic roots 
of (O have negative real parts and x = x* is locally asymptotically stable. □ 



From Theorems 13 . 21 and 13.31 it follows that the non-trivial equilibrium x = a;* of equation 
^ is locally asymptotically stable when 

We are going to show that as soon as condition does not hold, then the equilibrium can 
be destabilized. In the next section, we shall show that if condition PTI)l does not hold, then 
a Hopf bifurcation indeed occurs at x = x* . 



4 Hopf Bifurcations 

In this section we are going to show that the non-trivial equilibrium x = x* of equation ([T]) 
can be destabilized via Hopf bifurcations. The time delay r will be used as a bifurcation 
parameter. This result is obtained in Theorem 14. II 

Recall that the non-trivial equilibrium x = x* oi equation ([1]) exists if and only if (3o > 
(5 > 0. In the following, without loss of generality we assume that 

Again the results still hold when Tmin > 0, but the proof is easier to understand when 
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We look for purely imaginary roots of A(A). Of course, we assume that /3* < 0, otherwise 
a; = a;* is locally asymptotically stable. Let A ~ ioj, with a; G M, be a purely imaginary 
characteristic root of equation ([5]) . Then, t and uj satisfy the following system 



5 + (3*{1-2C{t,uo)) ^ 0, 
w + 2/3*S'(t,w) = 0, 



(13) 



where 

\ r 1 r 

C{t,uj) — — I coa{ujr)dr, S{t,ijj) — — sm{ujr)dr. 

T Jo T Jo 



First, one can see that iv = cannot be a solution of p3p . Otherwise 

S = (3* <0. 

Moreover, if w is a solution of system (|13p . then —ui is also a solution of Hence, we only 
look for positive solutions u. 

One can check that C(t, uj) and S{t, uj) are given, for r > and w > 0, by 

sin(wT) \ ^ 1 " cos(wt) 
C[T,UJ) ~ = K[ujt), j[t,uj) = , 

UJT UJT 

where the function K is defined by ([7])- Consequently, system (jl3p can be rewritten as 

K{^r) ^ (14) 



2/3* ' 
cos(a;T) — 1 1 
(wt)2 ~ 2/3*t' 



(15) 



Consider the sequence 
with 

In fact, one can check that 
Moreover, for all k eW , 



{xk]ken := {a; > ; x = tan(x)}, (16) 
Q — xq < xi < ■ ■ ■ < Xk < ■ ■ ■ ■ 

{xk}k&i = {2: > ; K'[x) = 0}. 

7r 

Xk e (fcTT, /CTT + -). 

Define two sequences {uk} and {ufc}, for fc e N, by 

Ufc := cos(x2fc+i) < 0, Wfc := cos(x2fc) > 0. 

Using the definition oi Xk, one can see that 

Uk = K{x2k+i) and Vk = K{x2k)- 

Thus, the sequence {ufcjfcgN is increasing with — 1 < < and the sequence {wfc}feeN is 

decreasing with vq = I and < Wfc < 1/2 for fc > 1 (see Figure [T]). 

Moreover, 

lim Uk ~ lim Vk ~ 0. 

k — >+oc k — > + oo 
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0.4 



8 10 12 14 16 18 



Figure 1: The graphe of K{x). 



Furthermore, one can check that, as soon as /3* < 0, 

S + (3* 



2p* 



<l = vo. 



Finally, define a function h, for x € [—1, 1/2), by 

2(1 -x) 



h{x) 



1 - 2x 



and set 



Hvq) = +00. 

We have the following results about the properties of the function h. 
Lemma 4.1. Suppose that 

00 -S 



h{uo) < n- 



f3o 



S + 13* ^0. 



(i) If 5 + 13* > 0, then there exists fc G N such that 



h{uk) < ^ < h{uk+i). 

(ii) If S + P* <0, then there exists k €N such that 

Hvk+i) < ^ < h{vk). 

Po 
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Proof. Since the function h is increasing on the interval [—1, 1/2), we can see that 



is equivalent to 



and 



is equivalent to 



h{uk) < < h{uk+i) 

Po 



6 + P* 

Uk < < Uk+l 



Hvk+i) < ^ < h{vk) 
Po 



Vk+l < < Vk. 



The lemma now follows. □ 
Proposition 4.1. (i) If 

f3 S 

h{uk) < n — < h{uk+i), /c G N, 

Po 

then system (T4\)- lil5\) has exactly 2(fc + 1) solutions {ti^i,uji^i), {Tk+i,i,LOk+i.i) o-nd 
(ri,2,wi,2); (Tfc+i,2,^fc+i.2) with 

0Ji,iTi^i e ((2/ - l)7r, X21-1), for I ^ 1, . . . ,k + I, 
0Ji,2n,2 G ix2i-i,2lir), forl^l,...,k + l 

and 

< n^l < ■ • • < Tfc+l_l < Tfe+1,2 < • • • < Tl_2. 

ft^) If 

n^^^h{uk), ken, 
Po 

then system \14-^- ll5]) has exactly 2fc + 1 solutions (Ti_i,a;i_i), {Tk+i,iT^k+i,i) and 

(n, 2,^1,2), (Tfc, 2,^^/0,2) with 

G ((2' - l)7r,a;2(-i), /or Z = 1, . . . , /c, 
'^i,2T;,2 e {X2i-i,2l7r), for 1 = 1,..., k, 

^k+l,lTk+l,l — X2k+1 

and 

< n.l < • • • < Tk+l,l < Tk,2 < ■ ■ ■ < Tl,2- 

(itl) If 

h{vk+i)<n^^-^ <h{vk), keW, 
Po 

then system has exactly 2fc + 1 solutions (ti_i, wi^i), (Tk+i,i,ijJk+i.i) and 

(n, 2,^1,2), (Tfc,2,Wfe,2) 

Wi,iTi,i e (7r/2,7r), 

Wi,iTi,i e (xi+i,(? + 2)7r), for I = 2,..., k + I, 

Wi,2Ti,2 e ((^ + l)7r, xi+i), for 1 = 1,..., k 
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and 

H If 



< < • • • < Tk+1,1 < Tk^2 < ■ ■ ■ < Tl,2- 

l3o-S 



f3o 



hivk), k G N*, 



then system ^T^-(T^ has exactly 2k solutions (ti^i, wi^i), (Tk^i,ujk^i) and (ti_2, ^1,2)7 

• • • , (Tfe,2,t^fc,2), with 

t^i,iTi,i G (7r/2,7r), 

Wi.iTi,! e {xi, {I + l)7r), forl^2,...,k, 
i^i.2Ti,2 e {{1 + 1)71, xi+i), /or / = 1, . . . , A; - 1, 

^k.2Tk.2 = X2k 



and 

(y)if 



< Tij < • • • < Tfe,! < rfe^2 < • • • < Tl,2- 



Po 



i/ien system ( [_?4[ )-|T5 )] /las a unique solution {ti,lui) such that ri > aKC? 

WiTi G (0, tt). 

Proof. We only prove (i) when fc = 0. The other cases can be deduced similarly. Assume 
that ^ 

h{uo) < — < h{ui). 
Po 



This is equivalent to 



S + f3* 



The function X is strictly negative and decreasing on (tt, xi) with K{y) G {uq, 0) (see Figure 
So the equation 

K{y) - 

has a unique solution yi on the interval (7r,a;i). Set 

(yi)^ 



'^1,1 



2/3*(cos(yi)-l) 



and — yi/ri^i. Then, (ti,i, cji^i) is a unique solution of system p4p - (|15p satisfying 
Wi,iTi,i G (7r,xi). 

Moreover, the function K is strictly negative and increasing on (a;i,27r) with K(ij) G 
(uo, 0), so the equation K{y) — [5 + (3*)/2(i* has a unique solution y2 on the interval (xi, 27r). 
Set 

^^'^ 2/3*(cos(y2)-l) 

and tJi 2 = y2/Ti,2- Then, (ti_2, 1^1,2) is a unique solution of system p4p - (fT5|) which satisfies 

Wl,2Tl,2 G {xi,2'k). 
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Furthermore, the function K is nonnegative on [0, tt] and 

Ui — K{x3) — min K{x). 

x>2tt 

Therefore, system (fn |) - ([TB1) has two solutions, (ti^i, wi^i) and (ti^2, '^i,2)- 
Finally, using the fact that 

cos(2/i) < cos(2/2), 

we obtain that 

This completes the proof. □ 

Lemma 14.11 and Proposition 14.11 give conditions for the existence of pairs of purely imag- 
inary roots of equation ([5|). In the next proposition, we study the properties of the purely 
imaginary roots of 

Proposition 4.2. Assume that there exists a Tc > such that equation (B^ has a pair of 
purely imaginary roots ±iujc for t — Tc with ujc > 0. // 

uJcTc 7^ Xk for all A; G N, 

where the sequence {xk}ken is defined by Mb]) , then zLiujc are simple roots such that 

dReiX) 



dT 
dRe{\) 



T — Tr 



> 0, if UJcTc & ix2k,X2k+l), 

< 0, if ujcTc e {x2k+i,x2k+2), fc e N. 

Proof. Assume that there exists a Tc > such that equation ([5|) has a pair of purely imaginary 
roots ±iujc for t = Tc with uJc > 0. Then, lUcTc satisfies system p^ - ([T5|) . 
Assume that 

oJcTc Xk for all k E N. 

Let us show that ztiWc are simple characteristic roots of ([5|). Using ([6]), one can see that 
iiwc are simple roots of (O if 

1 + 2(3* — {t,,uj,)^Q or (r„o.,)^0. 
ouj au) 

We will show that 

dC 

— (tcC^c) 7^ 0. 



A simple computation shows that 



dC . g{uJcTc 



where the function g is defined by 

g{x) = a;cos(a;) — sin(a;) for x > 0. 
One can check that g{x) = if and only if there exists a fco £ N such that x = Xk„ . Moreover, 
g{x) > if and only if x E {x2k+i,X2k+2), k eN. 
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This yields that 



and 



dC 



du 
du! 



{Tc,UJc)<0 if UJcTc e {x2k,X2k+l) 
{Tc,U)c)>0 if LOcTc e {x2k+l,X2k+2)- 



Hence, iiwc are simple characteristic roots of (0). 

Let A(r) — /i(r) + iuj{T) be a characteristic root of (O such that A(tc) = ±iu!c- By 
separating the real and imaginary parts, we obtain that 



Ai(r) + S + I3* 



2/3* 



cos(w(r)r)dr = 0, 



+ /" e-'^(^)''sin(a;(T)r)dr 



= 0. 



We denote by /^'(t) (respectively ij-''(t)) the first derivative of /i(T) (respectively w(t)) with 
respect to r. For r = Tc, we obtain that 



and 



r (99 
/x'(r,)[l + 2/3* — (r„c^,) 

9C 2/3* 

= 2/3* — (Tc,t^c)t^'(T-c) H (cOs(cJcT-c) -C(Tc,Wc)) 



r (99 

c^'(t,) 1 + 2/3* — (t„c^,) 
L aw 

dC 2/3* 
= -2/3* — (tc, Wc)/^'(rc) H (S{tc,uJc) -sin(wcrc)). 

OUJ Tr 



We consider two cases. First, assume that 



One can verify that 



1 + 2/3* — (re, We) = 0. 

OUJ 



1 + 2/3* — (r^c^e) = 2 + (5 + /3*)t,. 
OUJ 



Consequently, is equivalent to 



Then, it follows from equation (fT51) that 



5 + ^* 



dC I ^(tcCJc) - sin(wcTc) 

— (Tc^cj^i (Tc) = . 

Moreover, by using (fH]) and (fT5|) . wc have 

S'(tc, Wc) - sin(wcTc) _ 1 - ( cos(wcTc) + uJcTc sin(wcrc)) 



5 + (3* 
4/3* 



UJcT^ 



(17) 



(18) 



(19) 



(20) 
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Hence, (^01 implies that 



dC , ^ 



Since ^{tc,uJc) =/= 0, we have 



Furthermore, the sign of {tc) is the same as the sign of — ■|^(rc,Wc)- 



We now assume that 



1 + 2/3* — (^-^-)^0- 

OUJ 



Then, by using ([TT]) and (fTB]). we obtain that /^'(tc) satisfies 

dS , A 2 /„_9C n2 



2/3* 

Tr 



dC 

2/3* — (Tc,LJc)(S'(Tc,Wc) - sin(LJcTc)) 
+ (l + 2/3*|^(t„ We)) (cOs(cJerc) - C(Te, We)) 



Using the definitions of C and S, one can check that 
dC dS 

— {Tc,UJc){S{Tc,U!c) - sin(WcTc)) + — (Tc, We)(cOs(WcTc) - C{Tc,LUc)) = 0. 



Hence, 



M'(re)[(l + 2/3* — (re,We)) + (2/3* — (re, We) 

Notice that 



= ^^(cOs(Were) - C(Tc,We)). 



COs(WeTe) - C{Tc,UJc) _ gji^cTc) _ dC 



dS . 



Since 1 + 2/3* — (re, We) ^ 0, it follows that 



(l + 2/3*f(re,We))%(2r|^(re,We))%0 



dC , 



Consequently, fJ-'{Tc) ^ 0, and the sign of h'{tc) is the same as the sign of — ^(rc, We). 

In summary, we have obtained that, for t — Tc, equation ^ has a pair of simple purely 
imaginary roots itiwc such that 



dRe{\) 



dr 
dRe{X) 



dr 



> 0, 



<0, 



if UJcTc e (x2fe,X2fc+l), 
if Were e {x2k+l,X2k+2)- 



This completes the proof. 

Remark 1. If there exists a k €N* such that 



□ 



UJcTc = Xk, 
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then cither ±iLUc are not simple roots of ([5]) or 

dRe{X) 



0. 



dr 

Using a similar argument as in the proof of Proposition 14.11 we obtain that 



Thus, if 



5 + f3*' 

then HlOc are not simple roots of ([S]). If 
then 

dRe{X) ^ p 



dT 



r— Tc 



In the next theorem, we show that there exists a Hopf bifurcation at the non-trivial 
equilibrium a; = a;* of equation ([T]). 

Theorem 4.1. Assume that 

h{uo)<n^^-^ and 5 + 13*^0. 
Po 

Then a Hopf bifurcation occurs at x = x* for r = tq := min^^^^^^jj^.^ kefi '''c, where {tc, uJc) o,re 
solutions of |J^[ j- j75)) . defined in Proposition \4- 1\ When < t < tq, the equilibrium x = x* 
is locally asymptotically stable and it is unstable while tq < t < ti, where ti is the larger 
value of Tc such that lUcTc G {x2k,X2k+i), fc G N. 

Proof. We first check that a; = a;* is locally asymptotically stable when t £ [0,to). Notice 
that when t e [0,to), equation ^ does not have purely imaginary roots. Let r* > be 
small enough and fixed. Assume that, for r G (0,r*), equation ^ has a characteristic root 
A(t) = /Lt(r) + iuj{T) with /i(r) > 0. Separating the real and imaginary parts, we obtain 

Aifr) = -((5 + /?*) + ^ / e- coa(uj( T)r)dr 
T Jo 

and 

c^(r) = -^ re"^(")'-sin(w(T)r)dr. 



We deduce that, for r e (0, r*), 

|/i(T)| <\6 + l3*\^ 2/3* and \oj{t)\ < -2/?* 

Consequently, 

lim tii{t) — and lim tuj{t) — 0. 
Integrating by parts, we obtain 

2/x(t) = -{5 + /?*) + 2/3*e-^^(^)i^(™(r)). 
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Since /i(r) > 0, we have for r e (0, r*) that 

^{S + 13*) + 2/?*e-"'^(")X(ra;(r)) > 0. 
When T tends to zero, we obtain 

13* - S>0. 

However, /3* — S <Q. This is a contradiction. Therefore, for r G (0,r*), ii{t) < 0. Applying 
Rouche's Theorem [Sj p. 248], we obtain that all characteristic roots of (O have negative real 
parts when r € [0, min(rc)). Therefore, a; = a;* is locally asymptotically stable. 

Using Lemma 14.11 Propositions 14.11 and 14. 2|, we conclude to the existence of r/ . This 
concludes the proof. □ 

We illustrate the results of Theorem 14.11 in the next corollary. 

Corollary 4.1. Assume that the parameters S, (3q and n are given by Then there exists 
a unique value Tc > such that a Hopf bifurcation occurs at x = x* when t = Tc- When 
T < Tc, the equilibrium is locally asymptotically stable and becomes unstable when t > Tq. 
Moreover, when t = Tc, equation flj) has a periodic solution with a period close to 46 days 
(see Figure\^. The value of Tc is approximately given by 

Tc ~ 18 days. 

Proof. With the values given by Q , we obtain 

3 S 

~ 2.9153 > h(vi) ~ 2.3455. 

Po 

Hence, Proposition 14. II implies that the system (fTil) - P^ has a unique solution {tc,uJc) with 
Tc > and UcTc G (0,7r). From Theorem 14.11 we know that a Hopf bifurcation occurs at 
X = X* for T — Tc- The equilibrium is locally asymptotically stable when r < Tc and becomes 
unstable when r > Tc. Consequently, for r = Tc, equation ([T]) has a periodic solution with a 
period close to 2-k jtOc.. One can check that 

Tc ~ 18 days and lOc — 0.138. 

Computer simulations confirm our analysis (see Figure HJ. □ 

As mentioned earlier, the results in Theorem 14. II still hold when Tmin > 0. However, in 
this case, the computations in the proof of Theorem 14.11 are much more complicated. 



5 Discussion 

Hematological diseases have attracted a significant amount of modeling attention because a 
number of them are periodic in nature (Haurie et al. [8J). Some of these diseases involve 
only one blood cell type and are due to the destabilization of peripheral control mechanisms, 
e.g., periodic auto- immune hemolytic anemia (Belair et al. [2] and Mahaffy et al. |13j). 
Such periodic hematological diseases involve periods between two and four times the bone 
marrow production/maturation delay. Other periodic hematological diseases, such as cyclical 
neutropenia (Haurie et al. [S]), involve oscillations in all of the blood cells and very long period 
dynamics on the order of weeks to months (Fowler and Mackey [6 and Pujo-Menjouet et al. 
|14j ) and are thought to be due to a destabilization of the pluripotent stem cell compartment 
from which all types of mature blood cells are derived. 
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We have studied a scalar delay model that describes the dynamics of a pluripotent stem 
cell population involved in the blood production process in the bone marrow. The distributed 
delay describes the cell cycle duration. We established stability conditions for the model 
independent of the delay. We have also observed oscillations in the pluripotent stem cell 
population through Hopf bifurcations. With parameter values given in Mackey [9l[l0], our 
calculations indicate that the oscillatory pluripotent stem cell population involves a period 
of 45 days. 

It will be very interesting to study the dynamics of the two dimensional systems (Mackey 
[Ulin], Mackey et al. [12], Pujo-Menjouet et al. [H]) modehng the proliferating phase cells 
and resting phase cells with distributed delays. We leave this for future consideration. 



References 

[1] M. Adimy and F. Crauste, Global stability of a partial differential equation with dis- 
tributed delay due to cellular replication. Nonlinear Anal. 54, 8, 1469-1491 (2003). 

[2] J. Belair, M. C. Mackey and J. M. Mahaffy, Age-structured and two-delay models for 
erythropoiesis, Math. Biosci. 128, 317-346 (1995). 

[3] G. Bradford, B. Williams, R. Rossi and I. Bertoncello, Quiescence, cycling, and turnover 
in the primitive haematopoietic stem cell compartment, Exper. Hematol. 25, 445-453 
(1997). 

[4] F. J. Burns and I. F. Tannock, On the existence of a Gq phase in the cell cycle. Cell. 
Tissue Kmet. 19, 321-334 (1970). 

[5] J. Dieudonne, Foundations of Modern Analysis, Academic Press, New- York, 1960. 

[6] A. C. Fowler and M. C. Mackey, Relaxation oscillations in a class of delay differential 
equations, SI AM J. Appl. Math. 63, 299-323 (2002). 

[7] J. K. Hale and S. M. Verduyn Lunel, Introduction to Functional Differential Equations , 
Applied Mathematical Sciences 99, Springer- Verlag, New York, 1993. 

[8] C. Haurie, D. C. Dale and M. C. Mackey, Cyclical neutropenia and other periodic 
hematological diseases: A review of mechanisms and mathematical models. Blood 92, 
2629-2640 (1998). 

[9] M. C. Mackey, A unified hypothesis of the origin of aplastic anaemia and periodic 
hematopoiesis. Blood 51, 941-956 (1978). 

[10] M. C. Mackey, Dynamic hematological disorders of stem cell origin, in "Biophysical and 
Biochemical Information Transfer in Recognition," Eds. by J. G. Vassileva-Popova and 
E. V. Jensen, Plenum, New York, pp. 373-409, 1979. 

[11] M. C. Mackey, Cell kinetic status of haematopoietic stem cells, Cell Prolif. 34, 71-83 
(2001). 

[12] M. C. Mackey, C. Haurie and J. Belair, Cell replication and control, in "Nonlinear 
Dynamics in Physiology and Medicine", eds. by A, Beuter, L. Glass, M. C. Mackey and 
M. S. Titcombe, Springer, New York, pp. 233-269, 2003. 

[13] J. M. Mahaffy, J. Belair and M. C. Mackey, Hematopoietic model with moving boundary 
condition and state dependent delay, J. Theor. Biol. 190, 135-146 (1998). 



18 



M. Adimy, F. Crauste and S. Ruan 



A model of pluripotent stem cell dynamics 



[14] L. Pujo-Menjouet, S. Bernard and M. C. Mackey, Long period oscillations in a Go model 
of hematopoietic stem cells, J. Theoret. Biol, (submiited). 

[15] L. Pujo-Mcnjouct and M. C. Mackoy, Contribution to the study of periodic chronic 

myelogenous leukomia, Comptes Rendus Biologic (in press). 

[16] L. Sachs, The molecular control of hemopoiesis and leukomia, C. R. Acad. Sci. Paris 
316, 882-891 (1993). 



19 



